Chaos, Solitons & Fractals

As social beings, humans make decisions partly based on social interaction. Observing the behavior of others can lead to learning from and about them, potentially increasing trust and prompting trust-based behavioral changes. Observation-based decision making involves different neural structures. The orbitofrontal cortex (OFC) and lateral prefrontal cortex (LPFC) are known as neural structures mainly involved in processing emotional and cognitive decision values, respectively, while the anterior cingulate cortex (ACC) plays a pivotal role as a social hub, integrating the afferent expectancy signals from OFC and LPFC. This paper presents a neurocomputational model of the interplay between observational learning and trust, as well as their role in individual decision-making. Our model elucidates and predicts the emotional and rational behavioral changes of an individual influenced by observing the action-outcome association of an alleged expert. We have modeled the neurodynamics of three cortical structures (OFC, LPFC, and ACC) and their interactions, where the neural oscillatory properties, modeled with Dynamic Bayesian Probability, represent the observers attitude towards the expert and the decision options. As an example of an everyday behavioral situation related to climate change, we use the choice of transportation between home and work. The EEG-like simulation outputs from our model represent the presumed brain activity of an individual making such a choice, assuming the decision-maker is exposed to social information.

Voltage perturbations to a repetitively firing Hodgkin-Huxley (HH) model of neuronal spiking in the bistable regime with coexisting limit cycle and stable steady node can either lead to the spikes phase resetting or collapse to the stable steady state. The latter describes a non-firing hyperpolarized quiescent state of the neuron despite the presence of constant external current. Using asymptotic phase response curve (PRC), the impact of voltage perturbations on a repetitively firing HH model is studied here while it is diffusively coupled to another HH model under identical external stimulation. It is observed that the pre-perturbation state of synchronization and the coupling strength critically determine the PRC response of the perturbed HH dynamics. Higher coupling strengths of perfectly in-phase (anti-phase) synchronized HH models shrink (expand) the combinatorial space of perturbation strengths and the oscillation phases causing collapse to the quiescent state. This indicates reduced (enlarged) basin of attraction, viz. the null space, associated with the steady state in the HH phase space. The findings bear important implications to the spiking dynamics of diverse interneurons, as well as special cases of pyramidal neurons, coupled through electrical synapses via. gap junctions, and suggest the role of gap junction plasticity in tuning vulnerability to quiescent state in the presence of biological noise and spikelets.

Sex-related differences in the aging of the human brain were studied using large array of experimental data. The open archive CamCan was used as a source of data: the magnetic encephalograms, co-registered with magnetic resonance images of the head, were obtained for each of 434 subjects (ages 18-87 years, mean age 54.7 {+/-}18.4): 217 females (ages 18-87 years, mean age 54.5 {+/-}18.4) and 217 males (ages 18-84 years, mean age 54.8 {+/-}18.3). Recordings were split in 10-year age cohorts, each cohort consisted of equal number of men and women to calculate average intersex characteristics correctly. By massively solving the inverse problem, functional tomograms were calculated - the spatial distribution of elementary spectral components. Physiological noise was eliminated by joint analysis of MEG-based functional tomogram and magnetic resonance image for each subject. Then multichannel spectra were transformed into time series of the power of elementary current dipoles. Summary electric powers were calculated in six conventional frequency bands (1-4 Hz - delta; 4-8 Hz - theta; 8-13 Hz - alpha; 13-21 Hz - beta1; 21-30 Hz - beta2; 30-48 Hz - gamma), and sex differences in age-related changes were examined. It was found that in the youngest age cohort (18-29 years) the summary electrical power of the brain for males is 1.5 times greater than such power for females. For adults (30-69 years), male and female powers are approximately equal, while in older cohorts (70-87 years), male total brain power is greater. Age dependencies in various frequency bands are generally different for men and women, excluding higher frequencies 21-48 Hz. Basic conclusion can be made that after intersex averaging total electric power of the human brain is invariant through the lifespan from 18 to 87 years. The proposed method of joint MEG and MRI analysis can be used for further study of the sex-related details of brain sources in their connection with age changes.

Mass-conserved reaction-diffusion systems are used as mathematical models for various phenomena such as cell polarity. Numerical simulations of this system present transient dynamics in which multiple stripe patterns converge to spatially monotonic patterns. Previous studies indicated that the transient dynamics are driven by a mass conservation law and by variations in the amount of substance contained in each pattern, which we refer to as "pattern flux". However, it is challenging to mathematically investigate these pattern dynamics. In this study, we introduce a reaction-diffusion compartment model to investigate the pattern dynamics in view of the conservation law and the pattern flux. This model is defined on multiple intervals (compartments), and diffusive couplings are imposed on each boundary of the compartments. Corresponding to the transient dynamics in the original system, we consider the dynamics around stripe patterns in the compartment model. We derive ordinary differential equations describing the pattern dynamics of the compartment model and analyze the existence and stability of equilibria for the reduced ODE with respect to the boundary parameters. For a specific parameter setting, we obtained results consistent with previous studies. Moreover, we present that the stripe patterns in the compartment model are potentially stabilized by changing the parameter, which is not observed in the original system. We expect that the methodology developed in this paper is extendable to various directions, such as membrane-induced pattern control.

A novel mathematical framework to define the threshold of action potentials in excitable cells is presented. Unlike previously applied methods that rely on approximations or specific fixed-point bifurcations, the approach focuses on the geometry of membrane potential trajectories. Specifically, the focus is on the concavity changes during the upstroke of an electrical pulse. These changes in concavity form a curve of inflection points that defines a region in phase space crossed by all the action potentials in the system, and containing no non-action potential trajectories. Such region is called the excitability region and its size can be measured, thus providing a measure for the excitability of a dynamical system, and a way to compare the excitability between systems representing different biological phenotypes and stimulus conditions. The work transforms the traditionally vague physiological concept of excitability into a rigorous analytical description applicable across continuous, single compartment models of electrical excitability.

In this study, resource-constrained learning methods were developed as a model for the learning behavior of the fly brain, specifically the mushroom body. Recent research on the mushroom bodies of flies shows that unfamiliar odors activate certain output neurons (MBONs); however, these effects are rapidly suppressed upon repeated exposure to the same odor. Such MBON behaviors appear to reflect odor learning. We investigated how flies continue learning about odors throughout their lives despite their small brains. Researchers have suggested that learning about new odors can help flies forget existing memories. Therefore, we hypothesized that the main reason for continual learning is that it serves as a strategy for forgetting. To test the validity of this hypothesis, we designed three models using a kernel perceptron. This approach is suitable for estimating ongoing learning capacity within a budget. According to the results of computer simulations and theoretical analysis, the model demonstrated the importance of forgetting mechanisms for two reasons: first, to prepare for subsequent learning sessions, and second, to reduce the negative effects of deleting memories. Author summaryDrosophila mushroom body output neurons (MBONs) in the 3 compartment of the fruit fly brain are highly activated by novel odors, and their activation triggers alerting behavior. Interestingly, these specific neurons react only to unfamiliar odor information, suggesting they constantly undergo incremental learning of new odors. This study was aimed at constructing three incremental learning models of the MBON 3 neurons. Although there have been numerous studies on complex circuit designs to reproduce activation waveforms, herein we constructed a fundamental learning model based on a kernelized learning method. Since kernelized learning models interpret Hebbian learning as the addition or subtraction of kernel functions, the model is easy to analyze theoretically. Consequently, we conclude that the forgetting property observed in the MBON 3 neurons is essential for reducing error when learning occurs within a brain of limited capacity.

The brain is believed to process information efficiently in a different manner from deep learning-based artificial intelligence (AI). Brain-like next-generation AI is gaining attention owing to its potential to perform human-like, highly adaptive, robust, and power-efficient computation. To realize such AI, one crucial approach is the bottom-up implementation of the neuronal systems, capturing their electrophysiological characteristics in electronic circuits. However, this neuromorphic approach generally focuses on simplified neuronal models that do not refer to many biological findings. Developing closer-to-brain models is a natural direction that serve as a fundamental computing model for next-generation AI. One of the constraints of neuromorphic circuits is the bit resolution of synaptic efficacy memory, as the memory footprint scales with it precision. Although low-resolution synaptic efficacy is essential for minimizing memory circuit footprint and energy consumption, it generally leads to performance degradation in many tasks such as the spatio-temporal spike pattern detection. This study proposed a closer-to-brain learning rule that incorporates heterosynaptic plasticity (HP) induced by glutamate spillover. It is demonstrated that our model mitigates the performance degradation associated with low-bit resolution synaptic efficacy, achieving the pattern detection success rate with 3-bit resolution synaptic efficacy, which is comparable to 64-bit floating-point precision. Furthermore, the findings of the study indicate that HP based model accelerates the convergence of the synaptic effcacy and effectively potentiates the synapses relevant to the pattern detection while suppressing irrelevant ones, thereby promoting a bimodal distribution of synaptic efficacies. These findings may provide a basic framework for constructing an energy-efficient, brain-like next-generation AI that maintains high performance under hardware constraints.

BackgroundThe relevance of spontaneous activity has been unlocked thanks to recent large scale recordings that revealed, via Shared Variance Component Analysis (SVCA), the high-dimensional nature of the ongoing activity. A fundamental problem is how the dimension modifies when more neurons are included in the analysis. Contradictory results have been reported on this subject based on SVCA and Principal Component Analysis (PCA). New MethodWe investigate pro et contra of SVCA and PCA for the identification of reliable responses encoding underlying state variables. We focus on common features of the spectra of the reliable variances (RVs) and on their dimensionality. The analysis is demonstrated on previously published Ca2+ data from the visual and the dorsal cortex in head fixed mice during spontaneous behavior. ResultsRVs grow proportionally to the number N of neurons and show a power-law decay k- with the k-th SVC dimension over a range bounded by a maximal dimension kc, initially diverging as N 1/ and then saturating at sufficiently large N. The reliable dimensionality, estimated with different methodologies, also shows a clear saturation to an asymptotic value for large N. Furthermore, its value decreases when becomes larger, as demonstrated by employing experimental data as well as theoretical predictions. ConclusionWe have shown that SVCA is an extremely effective tool to extract reliable features from the neural signals, and that the exponent represents a biomarker able to reveal the level of correlation of the neurons as well as the dimensionality of the reliable space. HighlightsO_LIAdvantages and drawbacks of Shared Variance Component Analysis to extract reliable signals from neural data C_LIO_LIComparison of different methods to estimate reliable neural dimensionality associated to spontaneous activity C_LIO_LIAnalytical expressions of embedding dimensionality for power-law decaying reliable variances C_LIO_LIBounded growth of the dimensionality with the number of neurons C_LI

In previous studies by the author on binocular vision with the asymmetric eye (AE), which models a healthy human eye with misaligned optical components, the results were primarily presented in the Rodrigues vector (RV) framework and supported by simulations and 3D visualizations in GeoGebras dynamic geometry environment. In this paper, the novel geometric kinematics of the human eye, that is, the eye with misaligned optics, and simplified assumptions about the eye rotations (the eyes translational movements are disregarded), are developed within the framework of rigid-body rotations. The originality of the analysis lies in a precise geometric decomposition of a full rotation of the eyes posture into a torsion-free rotation (the geodesic part) and a torsional rotation (the non-geodesic extension of the geodesic part). This decomposition is extended to the corresponding decomposition of the angular velocity. A novel derivation of the eyes angular velocity from the RV formulation of the eye kinematics is proposed.

Computational models of brain structures impacted by Parkinsons disease are useful for exploring potential therapies. We use the Hindmarsh-Rose neuronal model to simulate synchronized activity in the subthalamic nucleus, capturing key features of the pathological rhythms observed in Parkinsons disease using a relatively small network of 100 neurons. Our model incorporates unidirectional excitatory chemical synapses whose strengths evolve according to a spike-timing-dependent plasticity (STDP) rule. To account for inputs from unmodelled neurons, both uniformly distributed white noise and Poisson noise were explored. White noise produced a single stable state of synchronized neuronal activity whereas Poisson noise resulted in two stable states, one synchronized and one desynchronized. We applied coordinated reset stimulation with a rapidly varying sequence (RVS CR) to examine its ability to reduce neuronal synchrony. The neuronal population was divided into subpopulations representing distinct physical sites of stimulation, as in deep brain stimulation therapy, and phase-shifted stimuli were delivered to each subpopulation in a random sequence. We explored how stimulation frequency and the number of stimulation sites affect the efficacy of RVS CR at desynchronizing the network. We demonstrate that RVS CR efficacy is sensitive to the depression-to-potentiation ratio in the STDP rule, which may be an important parameter to tune when reconciling simulations with experimental data. Numerical simulation of neuronal networks is constrained by computational resources when models demand large networks. This work proposes a model that demonstrates similar utility with a relatively small network, enabling researchers to study pathological neuronal activity and treatments more efficiently.

The phenomenon of splitting was originally observed in hamsters which, after prolonged exposure to constant light, exhibit two rest/wake cycles within a subjective day. Splitting is a consequence of the left and right suprachiasmatic nuclei (SCN) falling out of synchrony. While it is known that split activity is characterized by an antiphase relationship between the left and right SCN and between the core and shell within each hemisphere, the role of the commissural projections that connect the right and left SCN is not known. In the present study, we investigate the impact of the inter-hemispheric connections on the split and unsplit dynamics of a computational model of the bilateral SCN. Our model has 4 nodes corresponding to each right and left core and shell. We simulated our bilateral model under different lighting conditions and measured its period and the phase relationships among the 4 nodes. To further characterize the dynamics of the system, we performed a bifurcation analysis. We found that the bilateral model automatically splits unless entrained by bright light/dark cycles, or unless it has excitatory inter-hemispheric connections. This suggests that excitatory cross-connections may be important for freerunning behavior. We found that constant light of varying intensities transitions the model between split and unsplit activity only in very limited conditions, but the strength and polarity of the contralateral connections play a much greater role in this dynamical transition. These findings suggest that splitting may involve plasticity of the inter-hemispheric connections of the SCN.

Infectious disease surveillance systems in tropical countries show that respiratory disease incidence generally manifests as year-round activity with weak fluctuations and irregular seasonality. Previously, using a ten-year time series of influenza-like illness (ILI) collected from outpatient clinics in Ho Chi Minh City (HCMC), Vietnam, we found a combination of nonannual and annual signals driving these dynamics, but with unknown mechanisms. In this study, we use seven stochastic dynamical models incorporating humidity, temperature, and school term to investigate plausible mechanisms behind these annual and nonannual incidence trends. We use iterated filtering to fit the models and evaluate the models by comparing how well they replicate the combination of annual and nonannual signals. We find that a model including specific humidity, temperature, and school term best fits our observed data from HCMC and partially reproduces the irregular seasonality. The estimated effects from specific humidity and temperature on transmission are nonlinearly negative but weak. School dismissal is associated with decreased transmission, but also with low magnitude. Under these weak external drivers, we hypothesize that stochasticity makes a strong sub-annual cycle more likely to be observed in ILI disease dynamics. Our study shows a possible mechanism for respiratory disease dynamics in the tropics. When the external drivers are weak, the seasonality of respiratory disease dynamics is prone to the influence of stochasticity. Author SummaryAlthough the mechanisms driving seasonality of respiratory disease dynamics have been well-studied in temperate regions, they are unknown in the tropics. In this study, we used a 10-year influenza-like-illness (ILI) daily-reporting data set collected from outpatient clinics in Ho Chi Minh City (HCMC) in Vietnam to investigate the mechanisms associated with annual and nonannual ([~]215 days) periodic patterns in the data. By comparing seven mechanistic models against the data, we showed that the mechanism that best explains respiratory disease dynamics in HCMC is a stochastic susceptible-infected-recovered-susceptible (SIRS) model weakly driven by external drivers including specific humidity, temperature, and school term. The nonannual cycles duration is consistent with the inferred duration of immunity of the model. By showing the nonannual cycle as strong as in the data is only observed in stochastic model, we showed that the observed respiratory disease dynamics in HCMC is under the influence of stochasticity when external drivers are weak.

The flying fox bats roost in large colonies, suspended upside-down with minimal grip efforts from tree branches that are exposed to environmental disturbances. In this study, we investigate the oscillation dynamics of bats hanging from tree branches under natural conditions with wind. Bats modulate their grips to control the oscillation during wind disturbances and actively transform their postures. Using field observations, we analyze the angular deformation, speed, and phase of individual and collective bats swaying motions in response to environmental perturbations. We observed the mechanical coupling-based synchronization of collective bat oscillations on a tree branch. To rationalize this new phenomenon of bats synchronization behavior, we perform a table-top experiment of a physical model using active oscillators and passive systems. This work could inform the design of bio-inspired suspension systems and contribute to our understanding of animal balance and collective behavior in unsteady and complex environments.

Growing interest to describe the electrical behavior of glial cells, mainly astrocytes, in intact brain tissue poses more and more challenges to commonly accepted belief they only respond in a linear manner in uptake of the excess of extracellular potassium and maintenance of their network equipotentiality. Their highly conductive mutual interconnections via gap junction (GJ) connections introduce yet another class of nonlinear elements. As more studies report nonlinearities in membrane voltage Vm dependence of both, the membrane and junctional conductances, the need to formulate minimal dynamical models of their transient behavior is getting more acute. Since ODE models of coupled cells, even in simplest 1-d arrays, require simplified descriptions and small set of parameters, rare quantitative studies on glia makes the task even more difficult. This study attempts to qualify a self-coupled cell, or a glial cell coupled to fixed voltage as useful system for detecting the nature of instabilities and transitions coming from coupling. In a novel biophysical model of coupled astrocyte, we introduce nonlinear kinetics of deactivation for large junctional voltages for the first time. We found that N-shaped nonlinearities and corresponding fold structure in the vector field of isolated cell serves as a baseline on top of which coupling nonlinearities enrich the bifurcation picture. Numerical simulations of 1-d array of coupled astrocytes show that coupling increases the propensity of astrocytic Vm to bistability and front propagation. We believe that presented illustrations of possible effects of coupling nonlinearities will motivate neurobiologists to further explore their impact in disease. Significance statementTransient changes in membrane voltage of glial cells may produce significant transient voltage difference between directly coupled cells. Nonlinear steady-state conductance of their interconnection elements, the gap junctions, introduce nonlinear current profiles which are very difficult to measure and quantitate using the available methods due to marked permeability of the junctions and leakiness of glial membrane in general. We propose a minimal model of glial membrane extended with a self-coupled feedback loop, which under realistic simplifying assumptions could serve for qualitative analysis of the impact of coupling, on the stability of resting membrane voltage. Neuronal cells of the brain and spinal cord cannot exist and function without supportive and neuromodulatory functions of the diverse population of glial cells. This applies to virtually all physiological processes on cell level - from cell development, metabolic support, membrane signaling, slow molecular signal transduction, ion homeostasis, neurovascular coupling, myelination, to mention only a few, manifest neuro-glial interaction. Even though all glial cell types are interconnected, the most abundant ones, the astrocytes are massively interconnected by gap junctions to form ordered networks. Electrically, astrocytic networks display membrane voltage equipotentiality, which is considered system-wide resting state for given neuro-glial circuit or unit. With molecular and cellular substrates of glial connectivity being slowly elucidated, network science and dynamical modeling are slowly "invading" that area with many important issues left open. In this study using classical dynamical systems approaches we give indications how nonlinear intercellular coupling between astrocytes affects physiological resting state and its instabilities compared to isolated, uncoupled cell. We strongly believe the suggested minimal model could fill the gap in ODE modeling of neuro-glial circuits, within broadest scope of hypothesis-driven research in cell-level neuroscience.

The awake brain is known to display asynchronous (AS) states during periods of attention and arousal, but the responsiveness properties of such states remain unclear. Here, we investigate this question using computational models of spiking networks of excitatory and inhibitory neurons, mimicking recurrently-connected networks in layer 2/3 of the cerebral cortex. The networks can generate a continuum of AS states, but with different responsiveness characteristics. By using a mean-field model to infer the dynamic properties of the system, we find that there are two families of AS states, which we call "underdamped" (UD) and "overdamped" (OD). Responsiveness is maximised at the transition between OD and UD states, which identifies a "working point" that may present advantageous computational properties.

Alzheimers disease (AD) is characterized by the accumulation of amyloid-{beta} (A{beta}), yet the specific link between plaque burden and cognitive decline remains a subject of intense investigation. This paper presents a mathematical model that simulates the coupled dynamics of A{beta} monomers, soluble oligomers, and fibrillar species in the brain tissue. By modifying existing moment equations to include a dedicated conservation equation for A{beta} monomers, the model explores how various microscopic processes, such as primary nucleation, surface-catalyzed secondary nucleation, fibril elongation, and fragmentation, contribute to macroscopic disease progression. Central to this study is the concept of "accumulated neurotoxicity" as a surrogate marker of biological age, defined as the time-integrated concentration of soluble A{beta} oligomers. Unlike plaque burden, accumulated neurotoxicity cannot be reversed, and the harm it causes depends critically on the sequence of events that produced it. Numerical results demonstrate that while plaque burden and neurotoxicity both increase over time, their relationship is non-linear and highly sensitive to the efficiency of protein degradation machinery. Specifically, impaired degradation leads to a rapid advancement of biological age relative to calendar age. The model further identifies oligomer dissociation and fibril fragmentation as potential protective mechanisms that can counterintuitively reduce neurotoxic burden by diverting monomers away from the soluble oligomer pool. These findings provide a quantitative framework for understanding why individuals with similar plaque burdens may experience vastly different cognitive outcomes, underscoring the importance of targeting soluble oligomers early in therapeutic interventions.

The primary objective of system neuroscience is to understand the functional mapping and its causation in the dynamics of the brain network. Some experimental and methodological studies suggest that functional modularity and its hierarchical information processing in the brain network are crucial to understanding the functional role of task-specific or state-specific information flow in the brain. However, because most of the established techniques for detecting effective network structures in the neuroscience research field are strongly based on the "Granger causality" perspective, existing causal discovery methods specified for brain network analysis cannot identify the causal hierarchy in the modular network in the brain due to spurious correlation issues and indistinguishability of causal direction under the Gaussianity of observational noise in a linear system. To address the issues, we developed a causal discovery method for synchronous neural dynamics, called the Jacobian-informed linear non-Gaussian acyclic model, "j-VAR-LiNGAM", by incorporating the information of the Jacobian matrix determined from a phase-coupled oscillator model estimated from observed neural data into the VAR-LiNGAM algorithms. The method was validated by showing that it could extract causal ordering in both synthetic data and empirical neural observed data. Moreover, by analyzing the observed neural oscillatory signals obtained from mice and humans, we confirmed that our method identified causally hierarchical structures in the brain, which aligned with the neurophysiological interpretations. These findings suggested that our proposed method can reveal the neural basis of hierarchical information processing in the brain network.

BackgroundThis study examines a competition-based model (C-model) designed to capture the temporal dynamics of successive brain microstates derived from electroencephalography (EEG) recordings during eyes-open conditions. The analyzed data were obtained from a public repository comprising microstate sequences from 60 sessions of a single subject [1]. When applied to microstate dynamics, the C-model posits a stochastic competition among neural circuits underlying the expression of individual microstates. MethodsThe model is formulated at a conceptual level (computational level in Marrs framework) and employs a geometric distribution to account for the long right tail of microstate duration distributions, interpreted as the probability of "failure" of the currently active microstate to persist. To account for the short-lived left tail, the model incorporates a transient increase in the stability of the currently active network, or equivalently, a temporary decrease in the activation probability of competing microstates (refractory period). ResultsThe model provides a good fit to the microstate duration distributions across all 60 sessions. One third of sessions showed microstate identity sequential dependency with respect to the previous microstates. DiscussionThese results suggest that the C-model captures key aspects of microstate temporal structure. Moreover, because microstate probabilities can be modulated by psychophysiological conditions--including the influence of previously active networks--the model may serve as a building block for more comprehensive neurobiological frameworks of neural and behavioral dynamics. In such frameworks, microstate sequences could emerge from structured competition and flow among neural networks supporting microstate expression.

Sylvatic Yellow Fever (YF) is an infectious mosquito-borne disease with significant epidemiological relevance due to its widespread distribution and high lethality for human and non-human primates, particularly in tropical regions of the planet such as in Brazil. Identifying regions and periods of high environmental suitability for the occurrence of YF is essential for preventing or mitigating its burden, as it enables the efficient allocation of surveillance efforts, prevention, and implementation of control measures. Environmental modeling of YF occurrence has proven to be an effective approach toward this goal; however, its effectiveness strongly depends on the modeling framework's capabilities as well as the spatial and temporal precision of all associated data. We propose a fine-scale geospatial modeling of YF environmental suitability that is based on a generative machine-learning ensemble method built on a large set of high-resolution environmental covariates. First, we take the spatiotemporal statistical description of the environment of each of the 545 YF cases from 2019--2024 up to 30 m/monthly resolution at three buffer scales: 100 m, 500 m, and 1000 m ratios. Then, we perform a feature selection and train hundreds of One-Class Support Vector Machine submodels to form a robust ensemble model, whose predictions are projected to a 1x1 km resolution grid of Brazil under several metrics, exceeding seven million ensemble evaluations. The predictions ranked the Southern Brazil region with the highest mean suitability for YF, with a level of 0.64; Southeast comes next with 0.46, followed closely by Central-West region (0.44), North (0.39), and finally Northeast (0.28). The model exhibited high uncertainty for the North region, indicating that data collection efforts are much needed in this region. As for the environmental covariates, a feature analysis pointed out that Land use and cover accounts for the largest influence in the model output.

ContextGrid cells in the medial entorhinal cortex (MEC) of head-fixed mice exhibit ultraslow (<0.01 Hz) oscillations (USO) during walking in a 1D running wheel in darkness. It was proposed that these oscillations may have a connection with navigational behavior. ProblemThere is no clear link between the functional role of these oscillations and path integration, a fundamental navigation strategy used by animals to calculate their current position and orientation by continuously summing self-motion cues. HypothesisGiven the synaptic projections from MEC to the hippocampus, we hypothesized that ultraslow oscillations have a role in linking spatiotemporal memories acquired during navigation. MethodologyA realistic computational model of entorhinal-grid with ultraslow oscillations and hippocampal-place cells is proposed using synaptic plasticity between cell types, sustaining path integration of a rodent-like simulated animal. ResultsUltraslow oscillations induced persistent changes in the grid cell dynamics, represented as a positional drift of grid fields. Such drift resulted in position estimation errors but generated new grid-place cell associations when combined with synaptic plasticity. >DiscussionsUltraslow entorhinal oscillations were found to shape spatial memory through grid cell drifting, which could serve as a mechanism for flexibly accessing different spatial memories during navigation. HIGHLIGHTSO_LIPath integration dynamics hide ultraslow oscillations despite coexistence. C_LIO_LIUltraslow oscillations significantly degrade position estimation in path integration. C_LIO_LIGrid and place fields drift after the effect of ultraslow oscillations. C_LIO_LINew spatial memories were created as a result of the ultraslow oscillation drift. C_LIO_LIUltraslow oscillations enable dynamic access of different spatial memories C_LI